Universal contact and collective excitations of a strongly interacting Fermi gas 
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We study the relationship between Tan's contact parameter and the macroscopic dynamic prop- 
erties of an ultracold trapped gas, such as the frequencies of the collective oscillations and the 
propagation of sound in one-dimensional (ID) configurations. We find that the value of the contact, 
extracted from the most recent low-temperature measurements of the equation of state near uni- 
tarity, reproduces with accuracy the experimental values of the collective frequencies of the radial 
breathing mode at the lowest temperatures. The available experiment results for the ID sound 
velocities near unitarity are also investigated. 
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I. INTRODUCTION 

Significant experimental and theoretical work has been 
devoted in recent years to understand the universal prop- 
erties of interacting Fermi gases along the BEC-BCS 
crossover (for a review, see, for example, More re- 
cently, Tan has introduced a new concept for investigat- 
ing universality based on the so-called contact parameter, 
which relates the short-range features of these systems to 
their thermodynamic properties 0, Q . The universality 
of the Tan's relations has been proven in a series of ex- 
periments based on the measurement of the molecular 
fraction the momentum distribution, the RF spec- 
troscopic rate at high frequencies, the adiabatic sweep 
and virial theorems [5j] , the spin structure factor Q and 
the equation of state [7|. The temperature dependence 
of the contact parameter has been the object of recent 
theoretical [8] and experimental JOj work. 

In this paper we discuss the relationship between Tan's 
contact parameter and the frequencies of the collective 
oscillations of a harmonically trapped Fermi gas near 
unitarity. We also investigate the relationship with the 
sound velocity in highly elongated configurations. The 
study of the collective oscillations along the BEC-BCS 
crossover in terms of the contact parameter has been 
already addressed in [lo| where upper bounds to the 
collective frequencies were calculated using a sum rule 
approach. However, the sum rule method developed in 
this work significantly overestimates the hydrodynamic 
frequencies of trapped Fermi gases, being consequently 
ineffective for a useful quantitative comparison with ex- 
perimental data in the superfluid hydrodynamic regime 
of low temperatures. 

The present approach is based on a perturbative solu- 
tion of the hydrodynamic equations of superfluids near 
unitarity This allows for an exact analytic relation- 



ship between Tan'contact parameter and the deviations 
of the frequencies of the collective oscillations as well as 
of the one-dimensional (ID) sound velocity from their 
values at unitarity. The high precision achievable in the 
frequency measurements is in particular expected to pro- 
vide an alternative accurate determination of the contact 
parameter and to further confirm the universality of this 



physical quantity. 



II. CONTACT, EQUATION OF STATE AND 
HYDRODYNAMIC EQUATIONS 

We start from the following definition of the contact 
parameter, based on the so-called adiabatic sweep theo- 
rem 
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where E is the total energy of the system, I is the con- 
tact parameter, M is the atomic mass, and a is the s- 
wave scattering length. By using the local density ap- 
proximation (LDA), the total energy can be calculated 
a,s E — J d^r{e -\- nV^xt), where Vcxt is the trapping po- 
tential and e is the energy density of a uniform gas. The 
equilibrium profile, in the LDA, is determined by the 
equation 



where 



de(n, a) 
dn 
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(3) 



is the chemical potential of uniform matter, providing 
the equation of state of the gas, while fl is the chemical 
potential of the trapped system, fixed by the normaliza- 
tion condition. The derivative of the total energy with 
respect to 1 /a in Eq. ([1]) can then be conveniently written 
as 



(4) 



exploiting the link between the contact parameter and 
the equation of state. 

On the other hand the chemical potential of uniform 
matter /i(ri) is a crucial ingredient of the hydrodynamic 
equations of superfluids. At zero temperature, these 
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equations actually read [1/ 



at 



(5) 



= 0, 



so that an insightful question is to understand the link 
between the contact parameter and the behavior of the 
collective modes emerging from the solutions of Eq.(l5|). 
In the following we will discuss the problem near unitar- 
ity, where we can expand the chemical potential in the 
form 0, [Til 
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with ^ and C being two universal dimensionless param- 
eters accounting for the effects of the interactions. The 
first term in ^ exhibits the same density dependence 
as the ideal Fermi gas with the renormalization factor 
^. The second term (first-order correction in 1/a) is di- 
rectly related to Tan's parameter calculated at unitarity. 
Indeed, using Eqs.([l]|3]), one easily finds that the contact, 
for a harmonically trapped system, is given by 
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(7) 



where we have defined the Fermi wave vector kp — 
[37r^n(0)]^/'^ depending on the density in the center of 
the trap. 

The small-amplitude oscillations of the gas can be 
studied by solving the linearized hydrodynamic equations 
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where 5n is the amplitude of the density oscillations 
around the equilibrium value n. At unitarity oc 
m?/^, and the solutions of Eq.®, in the presence of har- 
monic trapping with axial symmetry, exhibit the analytic 
form 



^^W = ^ (4A2 + 5) ± v/l6A4 - 32A2 + 25 



(9) 



with A — ujz/uj^, holding for the lowest m — Q modes, 
and include the coupling between the monopole and 
quadrupole oscillations caused by the non spherical shape 
of the potential (here m is the z component of angular 
momentum carried by the excitation). Notice that, re- 
markably, result ([9]) does not depend on the parameter 
^ characterizing the equation of state at unitarity. This 
relation for the collective frequencies was actually first 
obtained in contexts different from the unitary Fermi 



gas, such as the ideal Bose gas above [13|, [14 1 
the ideal Fermi gas [Tsl in the hydrodynamic regime 



and 
In 



both cases the equation of state [^(ri, s) actually exhibits 



the same n^/^ dependence for fixed entropy per parti- 
cle s] and the hydrodynamic equations yield the same 
dispersion law ([9]) for the scaling solutions of coupled 
quadrupole monopole type in the presence of harmonic 
trapping. The prediction (|9]) for the collective frequen- 
cies was checked experimentally at unitarity providing a 
direct confirmation of the universality exhibited by the 
unitary Fermi gas [igI . \vf\ . The same oscillations were 
investigated out of unitarity along the BEC-BCS regime, 
confirming the predictions of theory [T^ [T^ and in par- 
ticular the fine details of the equation of state accounted 
for by quantum Monte Carlo simulations. 



III. COLLECTIVE FREQUENCIES AND 
SOUND VELOCITY SHIFTS NEAR UNITARITY 

In the following we calculate the deviations of the col- 
lective frequencies from the unitary value ([9]) holding for 
small values of the dimensionless parameter l/kpa. To 
this purpose we solve the hydrodynamic equations using 
a perturbative procedure, which is generally applicable 
to any equation of state having the form /x(n) = Hoiji) + 
/j,i(n), where /ii(n) ^ lji.vi{n) represents the first-order 
correction. The density profile n(r) = no(r) + n-i(i"), in- 
cluding the first-order correction, can be obtained start- 
ing from the expansion of ^{n) around the zero-order 
ground-state density profile no(r) and using the equilib- 
rium condition ^ in LDA. This gives 



m = [(5/i - /ii(no)] 



(10) 



where Sfi is the first-order correction to the chemical po- 
tential /I. Solving the linearized hydrodynamic equation 
(IH]) perturbatively, we find the following expression for 
the frequency shift: 



uj 2uj^M J d^rf* Sno 



(11) 



for the collective oscillations caused by the perturba- 
tion Hi{n) in the chemical potential. In Eg. (1111) /o — 
[diio{no)/dno]6no is the zero-order eigenfunction of ([8]), 
and ni is given by ([TU]). According to Eq. pTj) . one always 
has 6uj — for the surface modes satisfying the condition 
V^/o — 0, as expected, due to their independence of the 
equation of state. For compression modes one instead 
expects a correction due to the changes in the equation 
of state. In general, one can show that the first term 
in (llOp proportional to (5/1 gives no contribution to the 
frequency shift and will be consequently neglected in the 
following ^2^. 

Result pT|) is valid for both Bose and Fermi systems. 
In the case of weakly interacting Bose-Einstein condensed 
gas, it allows for the calculation of the frequency shifts 
caused by the Lee-Huang- Yang corrections in the equa- 
tion of state. In this case, the corresponding density cor- 
rection, neglecting the term proportional to 5/2, can be 
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FIG. 1: (Color online) Functions rj^ (solid blue) and 
(dashed red) relative to the higher and lower m — modes as 
a function of the deformation parameter A = oj^/ajx- 



written as rii = — 32(noa)'^/^/3-y/7r, and using Eg. ([TT|) . 
one finds the frequency shift of the compression mode as 
derived in [2l|. For the Fermi gas near unitarity, we in- 
stead employ the expansion of the equation of state © 
around the density profile no(r) calculated at unitarity. 
Ignoring also in this case the irrelevant term proportional 
to Sft, we find 



n{r) = no(r) + ni(r) = no(r) 



3/3 2/3 

— n, 
2a 



'0 



(r), 



(12) 



where a = fi^ ^{iir'^f I ^ / 2M , (3 = -2^\(37r2)i/3/(5Ma), 
and no = [(mo ~ Vl^xt)/^]'^^^, with /Jq being the chemical 
potential evaluated for a trapped system at unitarity. Af- 
ter some straightforward algebra, one finds the following 
expression for the frequency shift of the collective oscil- 
lations near unitarity: 



(13) 



One can check that this result is equivalent to the result 
of Eq.(24) and (25) in flj where a similar expansion was 
carried out near unitarity. 

The eigenfunctions for the m = modes EI) have the 
form fo^a + br\ 



cz^ , where 
Aio[4A2 + 3(a;/wx)2 - 10] 



b QMujI 
c ■i{uj/uji_f - 10 



(14) 



After some length but straightforward algebra, one fi- 
nally obtains the following result for the frequency shift: 
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U! 
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FIG. 2; (Color online) Frequency of the radial compression 
mode for an elongated Fermi gas in the units of radial fre- 
quency (cJx) [I3l- The deformation of the trap A = 0.038. 
The dashed red curve refers to the equation of state based on 
Monte Carlo simulations. Open and solid circles correspond 
to experiment measurements. The black straight line shows 
the slope of the frequency shifts around unitarity obtained 
from EQ. p5|) using ^ — 0.41 and ^ — 0.93 as extracted from 
the direct measurement of the equation of state 01 (corre- 
sponding to T/Nkp ~ 3.4). The green (light gray) straight 
line instead corresponds to the values ^ = 0.41 and = 0.74 
yielding X/Nk% ~ 2.7 (see text). 



with the index ± referring to the higher (-I-) and lower 
(— ) frequencies of ([9]). In the second equality of ([15]), 
we have used relation ([7]) for the contact parameter cal- 
culated for an harmonically trapped atomic cloud and 
we have expressed the Fermi momentum kp hi terms of 
the ideal Fermi gas wave vector k^p = (24A^)^/^a^J, with 
aho being the geometrical average of the harmonic os- 
cillator lengths. For the same total number of atoms, 
the density n(0) in the center of the trap for the unitary 
gas is ^^3/4 ^ijj^gg larger than for the ideal gas, yielding 
kp = kU-^l\ 

Eq. p^ represents the main result of the present paper. 
It relates Tan's contact, a central quantity for the univer- 
sality relations holding in interacting systems, with the 
low-energy macroscopic dynamics of the system, namely, 
the frequencies of the collective oscillations. These equa- 
tions can be used either to predict theoretically the fre- 
quency shifts, once the dimensionless parameters ^ and C 
or Tan's contact are known, or to determine experimen- 
tal constraints on the value of X. In FiglDwe plot ri± as 
a function of the deformation parameter A. For a spher- 
ical trap (A = 1) one obtains 774. = 1 for the monopole 
mode and rj- = 0, confirming, as already anticipated, 
that there is no frequency shift for the surface quadrupole 
mode. When A 9^ 1, the two modes are coupled. In the 
limit of both spherical traps and hi ghly elongated traps 
(A ^ 1) it reproduces the results of [ll| . 

In Fig[2]we show the prediction of Eq. p^ in the trap- 
ping conditions of the experiment of 17], using the values 
^ = 0.41, C = 0.93 extracted from the direct measure- 



4 



0.4 
0.36 
0.32 




-1 



-0.5 



0.5 



1/kOa 



FIG. 3: (Color online) Normalized ID sound velocity ci/vp vs 
the interaction parameter a in the unitary regime, where 
vf = hk%/M is the Fermi velocity for ideal Fermi gas and 
k% = fci?^^^** is the corresponding Fermi wave vector. Circles 
with error bars correspond to the experiment measurements 
[25j . The dashed red curve is a quadratic fitting [2^, and 
the black straight line is the slope of Sci/ci calculated from 
Ea. (|18|l using the values ^ — 0.41 and ^ = 0.93 extracted from 
0- 



ment of the equation of state Q carried out at the low- 
est temperatures and yielding the value I/Nkp ~ 3.4 
for the contact parameter. These values differ by a 
few percent from the most recent theoretical predictions 
based on Monte Carlo simulations at T = (see, for 
example, [2^ and references therein). The predicted 
slope Slu/uj ^ 0.11/(fc^a) (black line) of the collective 
frequencies of the radial breathing mode at unitarity 
turns out to be in very good agreement with experiments 
{Suj/uj - 0.12/(fc0,a) [23]). 

In order to appreciate the sensitivity of the slope to 
the choice of the values of ^ and in Figl^J we also show 
the predictions [green (light gray) line] for the frequency 
shifts using the values ^ = 0.41 and C = 0.74, yielding 
the smaller value I/Nkp ~ 2.7 for the contact. This 
value is closer to the measurement of the contact carried 
out in [sl and [oj at slightly higher values of temperature. 
The resulting slope Su/uj ^ 0.09/(fc^a) provides a worse 
description of the experimental data for the collective 
frequencies. 

The collective oscillations discussed above represent 
the discretized values of the usual sound waves described 
by hydrodynamics. It is actually useful to calculate also 
the changes of the sound velocity of a uniform sample 
near unitarity in terms of the dimensionless parameters 
^ and ^ or, equivalently. Tan's contact parameter. For 
bulk Fermi gases, one finds 5c/c — —(/{5£_kpa). 

In the case of cylindrical geometry with radial har- 
monic trapping the sound velocities can be also calcu- 
lated starting from the ID hydrodynamic expression (23 | 



Cl 




d r±n 



dfi{n) 



(17) 



should be evaluated in LDA. The above results hold for 
sound waves characterized by wavelengths significantly 
larger than the radial size of the gas. Carrying out a 
perturbative development around unitarity, similar to the 
one employed above for the calculation of the excitation 
frequencies, we obtain the result 
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where the the density profile along the radial direction 



where ci — {^/5y^^{hkp/M) is the ID sound velocity 
depending on the density in the center of the trap via 
kp. Notice that ci differs from the sound velocity cal- 
culated at the central value of the density by the factor 
■\/3/5. In FiglS] we show the experimental values of ci 
in the unitary regime measured in j25} together with the 
slope evaluated from Eq. ([T8)) using for £and ( the T ~ 
values obtained from the experiment [7|. A quadratic 
fit is applied to the experiment data (dashed red curve) 
[26j . The figure shows that the value of the slope at uni- 
tarity [black curve, Sci/ci ~ — 0.27/(fc^a)] overestimates 
the experimental linear changes in a visible way, sug- 
gesting that these experimental data were carried out at 
relatively higher temperatures, such that they cannot be 
accurately reproduced by employing the T = values of 
the contact parameter. Another source of disagreement 
might be due to the fact that the conditions of applica- 
bility of the ID expression (flT)) for the sound velocity are 
not fully satisfied in the experiment of [25} . 



IV. CONCLUSION 

In conclusion our analysis reveals consistency be- 
tween the experimental results for the contact, obtained 
through the measurement of the equation of state car- 
ried out at T ~ in [7], and the behavior of the collec- 
tive frequencies carried out in [l7j at the lowest temper- 
atures. It would be interesting to extend our analysis of 
the frequency shifts to finite temperature. The analysis 
could be simplified by the fact that the scaling modes of 
monopole and quadrupole types have a universal behav- 
ior at unitarity and their frequencies, calculated in the 
hydrodynamic regime, are independent of temperature. 
Furthermore they are not coupled to second sound. The 
proper calculation of the resulting slope at finite tem- 
perature and the explicit connection with Tan's contact 
parameter at finite temperature will be the object of a 
future work. 
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